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DATA FUSION OF STATIONARY ARRAY SENSOR AND SCANNING 

SENSOR MEASUREMENTS 

BACKGROUND OF THE INVENTION 
Field of the Invention 

[0001] This invention generally addresses the cross directional control of a 
process, such as a paper manufacturing process. The invention can improve the 
accuracy of a stationary array sensor employed for cross directional control by 
fusing the sensor output v^ith the output of a scanning gauge, or sensor, using a 
bank of filters such as Kalman filters, w^here each filter in the bank has a simple 
computational structure. 
Background Information 

[0002] Cross directional (CD) control of processes, such as paper 
manufacturing processes, is known. For example, U.S. Patent No. 4,903,528 
entitled "System and Process For Detecting Properties Of Traveling Sheets In 
Cross Direction", U.S. Patent No. 4,965,736 entitled "Cross-Directional Control 
Of Sheetmaking Systems" and U.S. Patent No. 5,121,332 entitled "Control System 
For Sheetmaking", the disclosures of which are all hereby incorporated by 
reference, are directed to cross directional control using a scanning sensor. A 
document entitled "Estimation of Cross-Directional Properties: Scanning vs. 
Stationary Sensors", Tyler, Matthew L. et al., AIChE Journal, Vol. 41, No. 4, 
April 1995, pages 846-854 also discusses cross-directional control of a process 
using a scanning sensor. The scanning sensor is used to perform process 
measurements in the cross direction of a moving sheet of paper. The 
measurements serve as feedback for control over some property, such as basis 



weight, moisture content or coating thickness, to render the property uniform 
across the moving sheet of paper. 

[0003] Because scanning sensors move back and forth across the paper sheet as 
it moves in a machine direction (MD), cross directional variations are not 
measured directly. Cross directional variations in a property to be measured can 
occur at a non-negligible rate relative to the speed of MD paper movement and the 
scan rate of the scanning sensor. Because the profile data associated with a 
scanning sensor is available only once per scan cycle, the scanning sensor's scan 
rate can be too slow to adequately address the dynamics of cross directional 
variations of the process property being controlled. 

[0004] Newer technologies have been proposed to provide faster sensing of the 
cross-directional properties of a product, such as the paper sheet in a paper 
manufacturing process, by using a stationary array sensor. A stationary array 
sensor includes a plurality of sensors located adjacent to one another in the cross 
direction, each sensor providing an approximately instantaneous measurement of a 
given property of the paper at a location across the paper's width. Although 
stationary array sensors avoid the need of a cross directional scan using motion of 
a single sensor back and forth over the moving paper sheet, the requirement that 
the stationary array sensor includes a plurality of individual sensors can render it 
quite expensive. For example, a stationary array sensor associated with a paper 
manufacturing process can require on the order of three hundred sensors to cover a 
sufficient width of the paper. This need for a large number of individual sensors 
in the stationary sensor array renders it cost prohibitive to include sensors having a 
high degree of precision. That is, to provide a stationary sensor array that can 
achieve measurements with an accuracy that approaches the accuracy of 
measurements achieved with a scanning sensor, expensive individual sensors must 
be used. 

[0005] It is also desirable to use stationary sensor arrays at locations in a 
manufacturing process which are unsuitable for scanning sensors. Again, in these 



circumstances, the stationary sensor array is typically configured with sensors that 
render the array cost competitive with scanning sensors used at other locations in 
the process. For example, stationary sensors are used in the earlier stages of a 
paper manufacturing process where the presence of flying debris, such as warm 
paper pulp, could jam the scanning mechanism of a scanning sensor. To achieve 
maximum benefit from the fast stationary array sensor measurement, such sensors 
would be used in close proximity to the actuators. Scanning sensors are commonly 
used in the end of the process, downstream from the actuator such that they 
measure properties of the finished product for quality control purposes. 
[0006] In addition to the use of less expensive, less accurate sensors in 
stationary sensor arrays, another factor which detracts from the quality of the 
measurements they provide is their susceptibility to drift. Although the sensors of 
stationary sensor arrays and scanning sensors can both experience drift, it is 
relatively easy to recalibrate a scanning sensor at least once during each scanning 
cycle. For example, the scanning sensor can be recalibrated during each cycle by 
moving it to a location off the paper being produced. Such a recalibration cannot 
be easily achieved with stationary sensor arrays, wherein each of the sensors is 
fixed in position. 

[0007] Thus, processes such as paper manufacturing processes and coating 
processes which involve cross-directional control, are known which use both 
stationary sensor arrays and scaxming sensors at various locations in the paper 
production process. It would be desirable to improve the accuracy of a stationary 
array sensor such that the quality of the measurement provided thereby is 
comparable to or exceeds that of a scanning sensor without rendering the stationary 
sensor array substantially more costly than a typical scanning sensor. 

SUMMARY OF THE INVENTION 
[0008] The present invention is directed to improving the accuracy with which 
a stationary array sensor provides cross directional measurements by periodically 
providing an offset compensation to the stationary array sensor using the output of 



a scanning sensor associated with the same manufacturing process. Exemplary 
embodiments correlate outputs from the stationary sensor array and the scanning 
array using a data reconciliation process. For example, a practical, real time data 
reconciliation of measurements from the scanning sensor and measurements from 
the stationary array sensor is achieved using a bank of Kalman filters to correlate 
outputs from the two sensors for each cross-directional measurement zone, 
wherein each filter possesses a relatively simple, computational structure. The 
Kalman filters can fuse the outputs from the stationary array sensor and the 
scanning sensor to track, and compensate, drift of the stationary array sensor. 
[0009] Generally speaking, the invention relates to a measurement system 
comprising: at least one stationary array of sensors at a first location to produce a 
first array of measurement outputs; at least one scanning sensor at a second 
location to produce a second array of measurement outputs; and means for 
synthesizing an array of measurement outputs by fusing (reconciling) the first and 
second arrays of measurement outputs. The invention is also directed at an 
associated method for fusing cross directional data measurements obtained from 
plural locations in a product manufacturing process. 

[0010] Exemplary embodiments compare and reconcile stationary array and 
scanning measurements so that the measurements can be correlated to the same 
spot on the manufactured material, such as paper. For example, measurements can 
be obtained which comprise time stamp information, cross direction coordinates, 
machine direction coordinates, and machine direction odometer or velocity 
information. 

[0011] The synthetic measurement can be obtained by computing a corrective 
offset (e.g., bias) updated by a recursive least mean square algorithm (e.g., update 
a bias model). For example, a filter, such as a bank of Kalman filters, can be used 
as the recursive least mean square algorithm to output data. The filter can be 
configured to compensate for different sensor inputs and bias errors. The filter 
can also compensate for the temporal variations in the biases of an array of 



stationary sensors. Data measurements obtained from stationary and scanning 
sensors can be compared by the filter, and an offset compensation for the drift of 
the stationary array sensor calculated. 

DESCRIPTION OF THE DRAWINGS 
[0012] The present invention will now be described by way of exemplary 
embodiments as illustrated in the following drawings: 
[0013] FIG. 1 is diagram of a paper production monitoring scheme and 
measurement setup; 

[0014] FIG. 2 is an overall system level diagram for the functions of the data 
fusion application; 

[0015] FIG. 3A illustrates one view of sensor measurement geometry; 
[0016] FIG. 3B illustrates another view of sensor measurement geometry; 
[0017] FIG. 4 illustrates a dynamic model for sensor bias; 
[0018] FIG. 5 illustrates a stationary sensor bias identification concept; and 
[0019] FIG, 6 illustrates moving-window least-squares processing to derive 
continuous CD paper parameter variations as a function of time, 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0020] The present invention relates to a measurement system, and associated 
method, which can be applied to any manufacturing process. For example, FIG. 1 
illustrates a measurement system for use with a paper production process wherein, 
according to exemplary embodiments of the present invention, data measurements 
obtained from plural locations using different types of sensors are fused. 
[0021] Referring to FIG. 1, a measurement system is illustrated for measuring a 
variable of at least one property of a product, such as a paper web, and including 
at least one stationary sensor and one scanning sensor. A stationary sensor shown 
in FIG. 1 can be a stationary sensor array 104 provided at a first location in the 
manufacturing process to produce a first array of measurement outputs. 
[0022] FIG. 1 also illustrates measuring the variable of the product and/or 
process with a scanning sensor such as scanning sensor 102. The scanning sensor 



102 is located at a second location in the manufacturing process to produce a 
second array of measurement outputs. The FIG. 1 measurement system includes 
means for synthesizing an array of measurement outputs by fusing the first and 
second arrays of measurement outputs using, for example, a synthesizing means 
represented in FIG. 1 as a processor 114. 

[0023] The processor 1 14 can fuse the outputs using a recursive least mean 
square algorithm implemented using a filter, such as a Kalman filter, to compute 
and/or update a corrective offset (e.e., update bias model). Although Kalman 
filters are known for fusing the measurements of different types of sensors, they 
have not been used to fuse the output of a stationary array sensor with a scanning 
sensor in manufacturing processes. The aforementioned Tyler et al paper mentions 
use of Kalman filtering when adding additional sensors to improve the estimation 
and control of a cross-directional property measured using a scaiming gauge. 
However, this document does not discuss using a scanning sensor to provide 
periodic off-set compensation for the measurement provided by a stationary array 
sensor, nor does it describe a Kalman filter which could be implemented in a cost- 
effective, practical fashion as a bank of multiple Kalman filters to achieve real 
time fusing of the outputs from these different types of sensors. 
[0024] FIG. 1 illustrates a web having an associated cross direction (CD) 103 
and a machine direction (MD) 105. An odometer 106 is used to measure where a 
point of the web L 101 has moved from a location associated with the stationary 
sensor array 104 to a location associated with the scanning sensor 102, and to 
provide MD odometer information. A transducer 107 can be used to output the 
position of the scanning sensor as it executes a CD traversal of the paper. Those 
skilled in the art will appreciate that the transducer can be located within the 
scanning sensor, however it can also be located at any other physical location 
where it is in operable communication with the scanning sensor. The odometer 
can be replaced by any device which can provide information representing a 



measure of movement including, but not limited to, velocity measurements using a 
web speed sensor and a dead reckoning algorithm. 

[0025] FIG. 2 shows an overall system-level diagram of data fusion application 
functions. The geometric registration application 207 and the model definition 
function 209, correspond to functions performed off-line. The middle two blocks 
depicting the geometric model 204 and the sensor measurement model 202 
describe parts of the model parameter database and correspond to static data that is 
not changing during the stationary operation of the manufacturing process. The 
model parameters define the internal dynamics and characteristics of the on-line 
Kalman filter application 210. Data from the scanner profiles 206 and the 
stationary array sensor profiles 208 is also input into the Kalman filter application 
210. The resultant output of this process is a synthetic measurement 212 that 
fuses, or merges, the measurements of the various sensor types found in the 
manufacturing process. 

[0026] The sensor measurement geometries that are available from an array of 
stationary sensors and a scanning sensor are depicted in FIGs. 3 A and 3B. The 
stationary array can include N sensors 301, spaced in the CD direction at a 
distance w 108. For example, in an exemplary paper manufacturing web process, 
each stationary sensor measures some paper characteristic, such as thickness, 
moisture content, coating thickness, or any other quality characteristic. These 
measurements are taken over a finite CD area, referred to as a data box 306 in 
FIG. 3 A. The stationary sensor measurements are compared with the 
corresponding measurements made by the scanning sensor. This comparison 
requires that all measurements be reconciled and synchronized so that 
measurements made by the two types of sensors are attributed to the same spot on 
the manufactured paper or the same point of whatever product is being measured. 
[0027] The measurements associated with each data box are an instantaneous 
snapshot in time comprising a large number of individual snapshots of the distinct 




-8- 

pixels that makeup the data box. These individual pixel measurements are 
averaged to form a single instantaneous measurement, 

[0028] Since the paper parameter being measured has no more than a linear 
spatial variation across the small data box in both CD and MD directions, the 
average instantaneous output of the stationary sensor will correspond to the 
parameter value located at the centroid of the data box. 

[0029] In FIG. 3 A, the paper can be divided into N zones 301 corresponding to 
the number of stationary sensors present in the system, with each zone being equal 
to the separation w 108, between the centerlines of the stationary sensors. At a 
downstream distance, L 101, from the line of stationary sensors, the scanning 
sensor traverses the paper, making a series of measurements of the same paper 
parameter of interest that the stationary sensor upstream of the stationary sensors 
had previously made. 

[0030] The scanning sensor output constitutes a series of discrete measurements 
of the paper parameter over small areas along a diagonal path 311 that is a function 
of the paper speed and the rate at which the scanning sensor is moving across the 
paper. As the scanning sensor crosses the different zones of the paper, a series of 
discrete measurements obtained in these regions is processed to derive the value 
corresponding to the midpoint of the zone. Purely random errors that occur in the 
array of stationary sensor outputs can be significantly attenuated. This is done by 
taking advantage of the characteristic of the paper parameters to vary in a smooth 
manner in CD, thereby allowing a moving-window least squares fit yielding a 
smoothed estimate of the paper parameter as a continuous function of CD position 
on the sheet, and as a continuous function of time. 

[0031] FIG. 3B shows another detailed view of the measurement data 
geometry. As the paper moves through the system the stationary array 
measurements 312 are taken. A fixed-point sensor measurement C-frame 
reference 308 is also taken. The sheet edge at the ends of the paper reels are 
shown by 314. There is also a region 304 where the scanner is moved off-sheet. 



One reason for moving the scanning sensor off-sheet can be for re-calibrating the 
sensor. Item 302 shows how the reel scanner measurement moves across the paper 
in a diagonal or zig-zag direction. 

[0032] In FIG. 3B, note that there are areas of missing measurements 300 that 
may occur. However these missing measurements will not prevent the present 
invention from working. A designer can empirically or experimentally determine 
what threshold level of missing measurements will resuU in reduced accuracy and 
compensate for these missing measurements. 

[0033] The invention permits direct comparisons to be made with the values 
measured by the appropriate stationary sensor at the center of that zone. A 
meaningful comparison would be made with respect to the value measured by the 
stationary sensor at an earlier point in time when the area of interest was at a 
distance L, upstream of the scanning sensor. The difference between the two 
sensor outputs may be used as the basis for continuous identification of the bias 
errors of the stationary sensors. 

[0034] Configuring a Kalman filter for fusing data measurements requires that 
an appropriate model for the temporal variations in each stationary sensor bias be 
defined. It is also useful to employ an approach that achieves maximum flexibility 
in representing possible variations of sensor bias. FIG. 4 shows a suitable model 
where two states are used to represent temporal variations in a parameter. For the 
special case where the state noise inputs r|l and t]2 are zero, the model reduces to 
a time-varying but non-stochastic function including a fixed offset plus a fixed 
drift. When the noise terms are activated, the model allows for random temporal 
variations in both the bias offset and bias drift. The selection of the spectral 
intensities of r| 1 and r|2 is based on the expected statistical variation in the bias 
resulting from temperature variations and other environmental and process factors 
affecting the variation. This selection affects the filter's ability to track real sensor 
bias variations without allowing it to become overly responsive. 
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[0035] A Kalman filter application executes on-line each time a new data 
arrives at the filter input, if there is a request for the filter output, or if there is a 
system timer request to run this module. 

[0036] The Kalman filter configuration can be based on a number of the explicit 
and implicit assumptions about the web process and the data measurements 
obtained. These assumptions lead to mathematical equations describing the 
problems to be resolved and the specific form of the algorithms required to solve 
the problems. These formalized assumptions will be presented as models of the 
system. The parameters of these models will be stored in an application database 
and used to define the tuning settings of these algorithms. The assumed models 
provide an abstraction layer interface between the process and the algorithm 
design. These models are but one exemplary embodiment of the present invention. 
Those skilled in the art will realize that other models can be used or the current 
models can be modified without departing from the spirit and scope of the present 
invention. 

[0037] The Kalman filter processes measurements of the differences between 
the quantities derived from the stationary sensors, and the corresponding quantities 
derived from the scanning sensor. The Kalman filter can be configured so that all 
of the information is fused (i.e., merged or integrated) in an optimal fashion. 
[0038] An effective Kalman Filter implementation includes the following 
elements. A model for dynamic state variations in the form of differential or 
difference equations is developed. Especially of interest in this exemplary 
embodiment are the errors found in the stationary sensor's bias compensation 
coefficients. In addition to the dynamic state variation model, a model for the 
random forcing functions is also developed. In this exemplary embodiment, the 
random inputs produce the random-like temporal variations of the stationary sensor 
biases . The other areas that must be modeled are the random errors that appear in 
any measurement equations that are used. In this exemplary embodiment, some 
examples of random error components are the outputs of the stationary sensor and 
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the scanning sensor, and measurement errors related to paper shrinkage and 
wandering as the paper moves from the Une of stationary sensors to the scanning 
sensor. 

[0039] An advantage of using a Kalman fiUer is that it can process all of the 
aforementioned data and model parameters in an organized and systematic way, 
thereby making it suitable for digital computer implementation. It also allows for 
convenient handling of non-uniform measurement sampling for each zone for the 
scanning sensor measurement. The following discussion summarizes the steps and 
mathematics that should be addressed in a Kalman filter implementation. 
[0040] Consider a system whose behavior is defined by the following set of 
discrete linear equations: 

where 

X = vector of states 

rj — vector of random (zero-mean) noise sequences 
^„ = state transition matrix from (n - 1)^^ to n^^ update points 
= noise distribution matrix 



[0041] For a given 0 and B, the state X will have a time variation determined 
by the particular noise sequence and initial condition, X^, which is generally 
taken to be a randomly distributed quantity. Since the noise sequence, has an 
infinite number of realizations, and the initial condition error can assume an 
infinite number of values, the system given by (1) has an infinite number of 
solutions. Because of this, attention is focused on the statistical behavior of 
Equation (1), rather than on specific solutions. 

[0042] A natural and useful way of characterizing the behavior of (1) is to 
compute the statistical parameters that define the bounds on the state vector, X. 
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The statistical bounds on the components of X are found by solving the covariance 
matrix equation associated with (1), which takes the recursive form: 

p p ^l+B Q (2) 



where P is the error covariance matrix of the state vector, X, defined explicitly by: 
and 

Pij=E(xpCj) 

[0043] in which E denotes the expectation operator. It is seen that the individual 
variances of the components of X are defined by the diagonal elements of P, with 
the joint expectations being defined by off-diagonal elements of P. The matrix Q 
in (2) is the covariance matrix of the driving noise vector, 77, defined by: 

Q - [qij 

in which 

[0044] Consider the case where the discrete process defined by (1) represents 
the true dynamic propagation characteristics associated with a given linear system. 
For this case, assume that a measurement is made at the nth measurement update 
time employing an external measuring device which allows a specific linear 
combination of the states to be directly monitored. A general way of stating this in 
mathematical terms is as follows: 
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n n ^ 



n 



(3) 



where 



= vector of measurements 
//„ = measurement matrix at nf^ measurement update time 
4 = measurement noise vector applicable to n^^ measurement 

and it is assumed that, in the general case, a number of independent measurements 

may become available simultaneously. 

[0045] The optimal utilization of information introduced through a series of 
measurements of the form given by (3), to estimate the state vector X in a 
sequential fashion, is the central problem addressed by Kalman estimation theory, 
and has the following solution. After each measurement (of a sequence of 
measurements), the estimate of the state, X, is refreshed by the two-step 
procedure: 



= optimal estimate of vector X just before the n^^ measurement is 

processed 



= optimal estimate of vector X hnmediately after the rf^ measurement is 

processed 

~ Kalman gain matrix at n^^ measurement update 



= ^X , 

n n n-\ 



(4) 



X = X +K\y -HX~] 



A 



(5) 



where 
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with being defined by 

K„ - P^HjiHf^Hj^R^y (6) 
in which 

~ apriori error covariance matrix of vector X 

= measurement noise error covariance matrix 
and the apriori error covariance matrix, P~ , is computed from (2) over the 
interval t , to t , 

[0046] After processing the ri^ measurement, the error covariance matrix of the 
state X is modified to reflect the benefit of incorporating new information 
introduced by the measurement as follows: 

P - {I-KH)P~ (7) 



where P^. is the aposteriori error covariance matrix. The form given by (7) is 
applicable when the Kalman filter is fully optimal; that is, when it is a full-state 
filter in which all components of X are fully accounted for in the mathematical 
model and, further, are re-estimated after each successive measurement is made 
available. 

[0047] Another issue is configuring a Kalman filter for the application of 
interest concerns the definition of an appropriate model for the temporal variations 
in each stationary sensor bias. An approach is desired that achieves maximum 
flexibility in representing possible variations in the sensor bias. The modeling of 
each stationary sensor bias is based upon a two state dynamic model that has been 
successfully used in other applications. The model is of sufficient generality to 
allow accurate tracking of slow temporal variations in the stationary sensors. This 
permits the high-frequency stationary sensor outputs to be compensated in real 
time to provide the information needed for effective process control. Because each 
Kalman filter consists of two states, the computational burden associated with 
cycling a large number of these filters is reasonable. 
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[0048] A model which has been successfully utilized in diverse applications 
uses two states to represent a temporal variation in a parameter. For the special 
case in which the state noise inputs (r|i, and r|2) are zero, the model reduces to a 
time-varying but non-stochastic function consisting of a fixed offset plus a fixed 
drift. When the noise terms are activated, the model allows for random temporal 
variations in both the bias offset and bias drift. The selection of the spectral 
intensities of rjj and 772 is based on the expected statistical variation in the bias 
resulting from variations in temperature and other environmental factors. This 
constitutes an important part of the filter specification, since it impacts the filter's 
ability to track real sensor bias variations without allowing it to become overly 
responsive, 

[0049] The state-space model for the stationary sensor bias variation is defined, 
with the aid of Figure 3, as follows 
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(8) 



[0050] Given the bias variation model defined by (8), the second major aspect 
of the Kalman filter that needs to be addressed is the measurement equation and 
measurement model. As discussed earlier, the output provided by the scanning 
sensor at the midpoint of each zone is used as the basis of a measurement if it is 
compared with the appropriate stationary sensor output at the same point on the 
paper. This is made precise by the following measurement equation 

ylsWkisyq'^is^L) (9) 



where 
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yf^(s) = measurement formed by comparing the output of the I^^ 

stationary sensor and the corresponding output of the scanning 
sensor at the midpoint of the zone 

s = distance traveled by the paper, as indicated by the odometer 
output 

L = distance between the array of stationary sensors and the scanning 
sensor 

q{(s) = output of the A:* stationary sensor at a point on the paper 
(s+L) = scanning sensor output at the mid-point of the kf zone at a 
point on the paper a distance L downstream 



[0051] Nominally, the two sensor outputs being compared are equal but, in 
reality, they will differ due to the stationary sensor bias, and the unavoidable 
random errors associated with the two sensor outputs. Therefore, the 
measurement error equation associated with the measurement equation defined by 
(9) is expressed in the time domain by 



3^/o=x,,(o+r-r (10) 



where 

x^^ = bias error in the I^^ stationary sensor at a time, t, corresponding to the 
odometer reading, s, that existed at the midpoint of the data box 
where the stationary sensor output occurred 

§^ = random measurement noise associated with the stationary sensor m 



random measurement noise associated with the scanning sensor 
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t — time corresponding to the paper distance s 

[0052] The variances assigned to the measurement errors, ^-^and are 
design parameters of the Kalman fiher that represent a tradeoff between the 
convergence characteristics of the parameter estimation process, and the filter's 
robustness in the presence of unmodeled (or unexpectedly large) measurement 
errors. The measurement errors account for the purely random errors in the 
sensor outputs, as well as errors introduced into the measurement from paper 
shrinkage and wander that can occur between the two sets of sensor outputs. The 
H matrix used in the Kalman filter update is defined by //=(1 0), 
[0053] The sensor measurement model will describe statistical characteristics of 
the measurement error in a compact form. It is assumed that the error is an output 
of a linear coloring filter driven by white Gaussian noise. The errors in the 
neighboring CD locations might be correlated with the cross covariance depending 
on the difference of the CD coordinates only and vanishing if this difference is 
large. In the initial stage, the coloring noise filter and the co variances will be set 
up empirically. At later stages, it might be possible to determine measurement 
models by analyzing on-line data and laboratory sensor calibration data. 
[0054] The process variation model describes the assumptions about the process 
change in time and CD coordinates. The process variation model will be shared 
among the different sensors (scanning and stationary) providing the measurements 
for the same process. An exemplary model is to assume that the measured paper 
property (e.g., weight) is random and independent in all points with different CD 
and MD coordinates. It is well-recognized that there is a correlation between the 
paper properties in the neighboring measurement locations. 
[0055] The geometric model defines relative position of the measurement 
datapoints obtained by different sensors as CD and MD. It can be assumed that 
the instantaneous measurement coordinates for the stationary array and the scanner 
are related through an affine transformation in CD coordinate and fixed MD offset. 
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A CD coordinate of each measurement can be computed through databox width 
sensor and a CD offset of this sensor. When computing the MD coordinates, 
additional machine (e.g. paper motion) speed and a time stamp for each 
measurement point can be taken into account. The relative CD offset of the 
sensors, the affine transformation defined by the paper sides wandering as it moves 
through the machine and the paper parameter variance (e.g. shrinkage) between the 
sensors. The MD distance between the sensors is not known accurately. 
Therefore these parameters need to be determined by processing data from an 
identification data set collected in a controlled experiment, such as a CD actuator 
bump test, to serve as a reference point. 

[0056] FIG. 5 shows one exemplary way of identifying biases in the array of 
stationary sensors 502. The outputs of stationary sensors 502 and a scanning 
sensor 516 are digitally pre-filtered to reduce noise before being utilized. A first- 
order lag filter would constitute a good pre-filter (518 for the stationary sensors 
and 516 for the scanning sensor). The same filter time constant should be used for 
all sensor outputs to avoid dynamic mismatch problems. 

[0057] An output from an odometer (106 of FIG. 1) controls the storage of the 
stationary sensor data in a moving-window memory 504. The storage interval 
should be selected to allow accurate interpolation of the discrete data in the later 
processing stages. This should not be too fine an interval, since this has a major 
impact on the required size of a storage device and the time required to refresh the 
memory in a moving- window fashion. 

[0058] The sensor outputs 520 from the moving window memory are processed 
522 by the respective Kalman filter in view of lagged stationary sensor outputs and 
measurement error. For example, the Kalman filter "1" update is applied to 
sensor "1" output. 

[0059] The moving-window memory is the means by which sensor outputs 
occurring at different points in time may be reconciled. The measurement of the 
paper characteristic obtained by the scanning sensor at the mid-point of a zone can 
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be compared to the corresponding output of the appropriate stationary sensor by 
simply going backward in the memory a distance L, using interpolation as required 
in order to arrive at synchronized measurements between the points. 
[0060] The moving-window memory stores the outputs of each stationary 
sensor at a specified increment of paper travel. The window size can be equal to, 
or greater than, the distance L between the array of fixed sensors and the line of 
travel of the scanning sensor. As each new set of sensor data is stored, the oldest 
set of data is eliminated from the memory and replaced by the new data. 
[0061] The scanning sensor position readout is used to control the selection of 
the Kalman filter to be updated based on the known zone width, w. 
[0062] The use of a bank of Kalman filters 501 to establish the gain and 
covariance matrices for each of the N stationary sensors allows the greatest 
flexibility in filter operation. This enables the system to deal with missing 
measurements, unequal update intervals for the various filters, failed sensors, etc. 
In FIG. 5 it is assumed that each filter is self-contained and can perform all the 
computations associated with the filtering process. It should be noted that the 
present invention suggests real-time operation, as each successive scanning sensor 
measurement becomes available. In practice, the scanning sensor may store and 
output data for a complete CD traversal rather than as a sequence of 
measurements. 

[0063] Additional memory for the collection of scanning sensor data is required 
and some additional logic allowing the measurements to be processed in the 
desired sequence needs to be added. The identification process defined above 
refi*eshes the bias compensation coefficients for the N stationary sensors at a rate 
equal to that of the scanning sensor as it traverses the paper in CD. This update 
rate is normally sufficient to allow tracking of the sensor bias temporal variations. 
Therefore, in a relatively short time the stationary sensor outputs will be 
accurately compensated for their bias errors. This leaves only the unavoidable 
purely random measurement error associated with each stationary sensor. 
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[0064] To address this error component, assume that the paper parameter of 
interest varies in a relatively smooth manner in the CD at any point along the 
direction of travel of the paper. Consequently, by assigning a polynomial function 
to the CD variation, and using the noisy, but bias-corrected outputs of the array of 
stationary sensors in least-squares fit, the effect of sensor noise can be attenuated. 
This is accomplished at the expense of a curve-fit error associated with the 
limitations of a given polynomial function to accurately represent the CD variation 
of the paper parameter of interest. 

[0065] A solution should balance the residual effect of sensor random 
measurement noise and curve-fit error is possible by defining a moving-window 
least-squares fit that Ihnits the CD window size over which the fitting function is 
used to represent the variations of paper parameters. In FIG. 6, the outputs of the 
stationary sensors 602(1), 602(2), 602(3) and 602(N) are input into their 
corresponding bias compensator 604(1), 604(2), 604(3) and 604(N). The output of 
the moving-window least-squares fit 606 is the function q(c,t) 608 which defines 
the paper parameter, q, as a continuous function of CD position, c, and time, t. 
[0066] Although the present invention has been shown and described with 
reference to exemplary embodiments, it will be understood by those skilled in the 
art that various other changes in the form and details may be made therein without 
departing from the spirit and scope of the invention. The invention can be used in 
any manufacturing process, in addition to the paper manufacturing process 
disclosed as an exemplary embodiment. 



